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For Coulomb blockade peaks in the linear conductance of a quantum dot, we study the correction to the 
spacing between the peaks due to dot-lead coupling. This coupling can affect measurements in which Coulomb 
blockade phenomena are used as a tool to probe the energy level structure of quantum dots. The electron- 
electron interactions in the quantum dot are described by the constant exchange and interaction (CEI) model 
while the single-particle properties are described by random matrix theory. We find analytic expressions for both 
the average and rms mesoscopic fluctuation of the correction. For a realistic value of the exchange interaction 
constant J s , the ensemble average correction to the peak spacing is two to three times smaller than that at J s = 0. 
As a function of J s , the average correction to the peak spacing for an even valley decreases monotonically, 
nonetheless staying positive. The rms fluctuation is of the same order as the average and weakly depends on J a . 
For a small fraction of quantum dots in the ensemble, therefore, the correction to the peak spacing for the even 
valley is negative. The correction to the spacing in the odd valleys is opposite in sign to that in the even valleys 
and equal in magnitude. These results are robust with respect to the choice of the random matrix ensemble or 
change in parameters such as charging energy, mean level spacing, or temperature. 

PACS numbers: 73.23.Hk, 73.40.Gk, 73.23.-b 



I. INTRODUCTION 

Progress in nanoscale fabrication techniques has made pos- 
sible not only the creation of more sophisticated devices but 
also greater control over their properties. Electron systems 
confined to small regions - quantum dots (QD) - and espe- 
cially their transport properties have been studied extensively 
for the last decade. 1,2 One of the most popular devices is a lat- 
eral quantum dot, formed by depleting the two-dimensional 
electron gas (2DEG) at the interface of a semiconductor het- 
erostructure. By appropriately tuning negative potentials on 
the metal surface gates, one can control the QD size, the num- 
ber of electrons n it contains, as well as the tunnel barrier 
heights between the QD and the large 2DEG regions, which 
act as leads. Applying bias voltage V between these leads 
allows one to study transport properties of a single electron 
transistor (SET), Fig.GJa).! 

We study properties of the conductance G through a QD in 
the linear response regime. We assume that the dot is weakly 
coupled to the leads: Gl,r -C e 2 //i, where Gl,r are the con- 
ductances of the dot-lead tunnel barriers, e > is the elemen- 
tary charge, and h is Planck's constant. 

To tunnel onto the quantum dot, an electron in the left lead 
has to overcome a charging energy Ec — e 2 /2C, where C is 
the capacitance of the QD, a phenomenon called the Coulomb 
blockade. However, if we apply voltage V g to an additional 
back-gate capacitively coupled to the QD, see Fig. [Da), the 
Coulomb blockade can be lifted. Indeed, by changing V g one 
can change the electrostatics so that energies of the quantum 
dot with n and n + 1 electrons become equal, and so an elec- 
tron can freely jump from the left lead onto the QD and then 
out to the right lead. Thus, a current event has occurred, and a 
peak in the conductance corresponding to that back-gate volt- 
age, V g n+ i, is observed. By sweeping the back-gate voltage, 
a series of peaks is observed, Fig.^b). 

In this paper we calculate the correction to the spacing be- 
tween Coulomb blockade peaks due to finite dot-lead tunnel 



coupling. In recent years, low-temperature Coulomb block- 
ade experiments have been repeatedly used to probe the en- 
ergy level structure of quantum dots. 1 - 2 - 5 The dot-lead tunnel 
coupling discussed here may influence such a measurement 
- the presence of leads may change what one sees - and so 
an understanding of coupling effects is needed. One dramatic 
consequence is the Kondo effect in quantum dots. 3 - 4 Here we 
assume that T 3> Tk, where Tk is the Kondo temperature, 
and, therefore, do not consider Kondo physics, focusing in- 
stead on less dramatic effects that, however, survive to higher 
temperature. 

We study an ensemble of chaotic ballistic (or chaotic dis- 
ordered) quantum dots with large dimensionless conductance, 
9^1. The dimensionless conductance is defined as the ratio 
of the Thouless energy Et to the mean level spacing A: g = 
Et/A> 5 For isolated quantum dots with large dimensionless 
conductance, the distribution of g energy levels {Ek} near the 
Fermi level and the corresponding wave functions {^(r)} 
can be approximated by random matrix theory (RMT).5i&i As 
will be evident from what follows, the leading contribution 
to the results obtained here comes from £ = 2Ec / A energy 
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FIG. 1: (a) Scheme of the Coulomb blockade setup; (b) Oscilla- 
tions of the SET linear conductance G as the back-gate voltage V g 
is changed. "Even" ("odd") corresponds to an even (odd) number of 
electrons in the valley. Arrows show average shifts in the positions 
of the peaks' maxima due to the finite dot-lead couplings. 
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levels near the Fermi level; thus, if Ec ^ £t the statistics 
of these £ levels can be described by RMT. We furthermore 
neglect the spin-orbit interaction and, therefore, consider only 
the Gaussian orthogonal (GOE) and Gaussian unitary (GUE) 
ensembles of random matrices. 

The microscopic theory of electron-electron interactions in 
a quantum dot with large dimensionless conductance brings 
about a remarkable result. 8,9 To leading order, the interac- 
tion Hamiltonian depends only on the squares of the fol- 
lowing two operators: (i) the total electron number operator 
h = J^cl Cfe CT where {cka} are the electron annihilation 
operators and a labels spin, and (ii) the total spin operator 
S = | J2 c Li I & \ G ^I c kcr 2 where {ct,} are the Pauli ma- 
trices. The leading-order part of the Hamiltonian reads as^ 
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c n 



(1) 



where Ec is the redefined value of the charging energy 8 ^ and 
J s > is the exchange interaction constant. Higher-order cor- 
rections are of order A / g&iLIiLlL The coupling constants in Q 
are invariant with respect to different realizations of the quan- 
tum dot potential. This "universal" Hamiltonian is also invari- 
ant under arbitrary rotation of the basis and, therefore, com- 
patible with RMT. In principal, the operator of interactions in 
the Cooper channel can appear in the "universal" Hamiltonian 
for the GOE case. However, if the quantum dot is in the nor- 
mal state at T = 0, then the corresponding coupling constant 
is positive and is renormalized to a very small valueAi^ The 
"universal" part of the Hamiltonian given by Eq. Q is called 
the constant exchange and interaction (CEI) model>I2iAi 

The total Hamiltonian of the quantum dot in the g — > oo 
limit thus has two parts, the single-particle RMT Hamiltonian 
and the CEI model describing the interactions. The capaci- 
tive QD-backgate coupling generates an additional term that 
is linear in the number of electrons: 
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where J\f = C g V g /e is the dimensionless back-gate voltage 
and C g is the QD-backgate capacitance. 

The CEI model contains an additional exchange interac- 
tion term as compared to the conventional constant interac- 
tion model (CI model). 113 Exchange is important as J s is of 
the same order as A, the mean single-particle level spacing. 
Indeed, in the realistic case of a 2DEG in a GaAs/AlGaAs 
heterostructure with gas parameter r s = 1.5, the static ran- 
dom phase approximation gives J s « 0.31 A. 14 Therefore, as 
we sweep the back-gate voltage, adding electrons to the quan- 
tum dot, the conventional up-down filling sequence may be 
violatediA£ii& Indeed, energy level spacings do fluctuate: If for 
an even number of electrons n in the QD the corresponding 
spacing, £a+i — en., is less than 2J S , then it becomes en- 
ergetically favorable to promote an electron to the next or- 
bital instead of putting it in the same one; thus, a triplet state 
(S = 1) is formed. Higher spin states are possible as well. For 
r s = 1.5 the probability of forming a higher spin ground state 
is P^S > 0) » 0.26 and P 2 (S > 0) » 0.19 for the GOE 
and GUE, respectively. The lower the electron density in the 



QD, the larger r s and, consequently, the larger the exchange 
interaction constant J s . 

The back-gate voltage corresponding to the conductance 
peak maximum M n _i is found by equating the energy for 

n—l electrons in the dot with that for n electrons: 17 



(3) 



As we are interested in the effect of dot-lead coupling on these 
peak positions, it is natural to expand the energies perturba- 
tively in this coupling: E = E^ + E^ + . . . . One possible 
virtual process contributing to E^ is shown in Fig.|2] Elec- 
tron occupations of the QD "to the left" and "to the right" of 
the conductance peak [see Fig.[3b)] are different; hence, the 
corrections to the energies are different. Therefore, the 
position of the peak maximum acquires corrections as well, 
J\f = A/" (0) + N {2) + . . . , as does the spacing between two 
adjacent peaks. 

This physical scenario has been considered by Kaminski 
and Glazman with the interactions treated in the CI model, 
i.e. neglecting exchange. 18 The ensemble-averaged change 
in the spacing and its rms due to mesoscopic fluctuations 
were calculated. On average, "even" spacings (that is, spac- 
ings corresponding to an even number of electrons in the val- 
ley) increase, while "odd" spacings decrease (by the same 
amount) ^ 
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where U is the dimensionless spacing normalized by 2Ec and 
9l,r — Gl.r/ [2e 2 1 h) are the dimensionless dot-lead conduc- 
tances. 

In this paper we calculate the same quantities, but with the 
electron-electron interactions in the QD described by the more 
realistic CEI model. We find that the average change in the 
spacing between conductance peaks is significantly less than 
that predicted by the CI model. However, the fluctuations are 
of the same order. In contrast to the CI result, 18 for large 
enough J s , we find that "even" spacings do not necessarily in- 
crease (likewise, "odd" spacings do not necessarily decrease). 

The paper is organized as follows. In Sec.|II|we write down 
the total Hamiltonian of the system and find the condition for 
the tunneling Hamiltonian to be considered as a perturbation. 
In Sec.|in]we describe the approach and make symmetry re- 
marks. In Sec. IIVI we perform a detailed calculation of the 
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FIG. 2: One example of the virtual processes contributing to E' 2 ' . 
This virtual process corresponds to an electron tunneling out of the 
quantum dot into the left lead and then tunneling back into the same 
level in the QD. 
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correction to the spacing between Coulomb blockade peaks 
for the i — > 1 — > | spin sequence. In Sec.|V]we find the 
ensemble-averaged correction to the peak spacing. The rms 
of the fluctuations of the correction to the peak spacing is cal- 
culated in Sec. lVIl In Sec. lVIII we summarize our findings and 
discuss their relevance to the available experimental data»A2i2S 



II. THE HAMILTONIAN 

The Hamiltonian of the system in Fig.^a) consists of the 
QD Hamiltonian [Eq. 10], the leads Hamiltonian, and the tun- 
neling Hamiltonian accounting for the dot-lead coupling: 



H = Hdot + H\ ca ds 



Hf. 



The leads Hamiltonian can be written as follows: 



loads 



pa 



e q C qcr C 



(5) 



(6) 



where {e p } and {e q } are the one-particle energies in the left 
and right leads, respectively, measured with respect to the 
chemical potential (see Fig.|2j. We assume that the leads are 
large; therefore we (i) neglect electron-electron interactions in 
the leads and (ii) assume a continuum of states in each lead. 
The tunneling Hamiltonian is 21 

Htun = ^(ifep4 CT Cp (T +h.C.)+^(t/cqcj £( T C q^+ h - C -), ( 7 ) 
kpa fcqo - 

where {tk P } and {ifc q } are the tunneling matrix elements. 

We assume that T<A and, therefore, neglect excited states 
of the QD concentrating on ground state properties only. We 
also assume that the QD is weakly coupled to the leads, treat- 
ing the tunneling Hamiltonian as a perturbation. Corrections 
to the position of the peak maximum can be expressed in terms 
of corrections to the ground state energies of the QD via Eq. 
0. The perturbation series for these corrections contains only 
even powers as H tua is off-diagonal in the eigenbasis of Hq. 
The 277i th correction to the position of the peak is roughly 
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Thus, finite-order perturbation theory is applicable if 
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To loosen this restriction one should deploy a renormalization 
group technique which, however, is beyond the scope of this 



paper. 
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III. PLAN OF THE CALCULATION 

As the exchange interaction constant J s becomes larger, 
more values of the QD spin S become accessible. The struc- 
ture of the corrections to the ground state energies depends 



on the total QD spin S, and this structure becomes very com- 
plicated for large values of the spin S. Fortunately, for the 
realistic case r s — 1.5, the probability of spin values higher 
than i in an "odd" valley is small: P2(S > i) ~ 0.01 for the 
GUE. Hence, we can safely assume that in the "odd" valley 
the spin is always equal to i. In the "even" valley, one has to 
allow both 5 = and 5=1 states. 

The structure of the expression for the spacing between 
peaks depends on the allowed spin sequences. For an "even" 
valley there are only two possible spin sequences: 







and 



1 



(10) 



where the number in the middle is the spin in the "even" val- 
ley, while the numbers to the left and right are spin values in 
the adjacent valleys, Fig.^b). For an "odd" valley there are 
four possibilities: 
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To obtain correct expressions for the average spacing between 
peaks, one should weight these sequences with the appropriate 
probability of occurrence. 

Before proceeding with the calculations, we note several 
general properties. First, ensemble-averaged corrections to 
the "odd" and "even" spacings are of the same magnitude and 
opposite sign, Fig.^b). Second, the mesoscopic fluctuations 
of both corrections are equal. Indeed, the shift in position of 
an "even-odd" (n — > n+1) peak maximum, Fig.^b), is de- 
termined by the interplay between the — > | and 1 — ► | spin 
sequences. Likewise, the shift of the "odd-even" (77 — 1 — > 77) 
peak is determined by the | — > and 5 — ► 1 spin sequences. 
Now if we sweep the back-gate voltage in the opposite direc- 
tion and write the same peak as 77 — > 77—1, then the corre- 
sponding spin sequences are exactly the same as they were in 
the first case: — > | and 1 — ► |. From this symmetry argu- 
ment one can conclude that (i) the ensemble-averaged shifts 
of the "even-odd" and "odd-even" peaks are of the same mag- 
nitude and in the opposite directions 18 and (ii) the mesoscopic 
fluctuations of both shifts are equal. 

Thus, to simplify the calculations we study only the "even" 
spacing case. This corresponds to the two spin sequences 
given in Eq. fllOb . First, we calculate corrections to the spacing 
between peaks for both spin sequences. A complete calcula- 
tion for the doublet-triplet-doublet spin sequence is in the next 
section. Second, we elaborate on how to put these spacings to- 
gether in the final expression for an "even" spacing. Finally, 
we calculate GOE and GUE ensemble-averaged corrections 
to the spacing and the rms fluctuations. 



IV. DOUBLET-TRIPLET-DOUBLET SPIN SEQUENCE: 
CALCULATION OF THE SPACING BETWEEN PEAKS 



Let us find the correction to the spacing between peaks for 
a doublet- triplet-doublet spin sequence. The corresponding 
electron occupation of the quantum dot in three consecutive 
valleys with n— 1, n, and n+1 electrons is shown in Fig. [3] 
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For the isolated QD, the position of the n—1 
tance peak maximum is determined by 



n conduc- Thus the spacing between peaks in zeroth order is 
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where fl=i and E^' s= , are the ground state energies of 

the dot Hamiltonian [Eq. (|2)]. The corrections due to dot-lead 
tunneling are different for the doublet and triplet states. The 
resultant shift in peak position is given byi& 



2Er 



Note that for the second-order correction to the position, the 
ground state energies are taken at the gate voltage obtained in 
the zeroth-order calculation, Eq. dl 2i . 

Analogous equations hold for the n — > n + 1 conductance 
peak. The spacing between these two conductance peaks is 
then defined as 

Z4, s =i =K+d 1 - \) -> !)• ( 14 ) 



A. Zeroth Order: Isolated Quantum Dot 



For the doublet-triplet n — 1 — > n sequence, Eq. (I12> gives 



1 

"~2 



2Er 



4 2 2 



(15) 

where the last temperature-dependent term is the entropic cor- 
rection to the condition of equal energies. 24 For the position 
of the n — » n+1 peak maximum we obtain 
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where j = J s /A and a; = — £f) /A (see Fig. [3j. Sim- 

ilarly, for the doublet-singlet-doublet spin sequence the spac- 
ing is 
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(13) Note that in both cases depends only on the spacing x. 



B. Second Order: Contribution From Virtual Processes 

Let us consider in detail the second-order correction to the 
ground state energy of the triplet for subsequent use in Jl 3i : 
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where the sum is over all possible virtual states. E^ and 
are the eigenvalues and eigenvectors of H^ ot , Eq. (|2}. 

Different terms in Eq. (1191 have a different structure de- 
pending on the type of virtual state involved; six possibilities 
are shown in Fig. [4] To take into account all virtual processes, 
we sum over all energy levels in the QD and integrate over 
states in each lead. To simplify the calculation even further, 
we assume (just for a moment) that T = so that the Fermi 
distribution in the leads is a step function. Later, we will see 
how T reappears as a lower cutoff within a logarithm. 

Following the order of terms in Fig. |4j the second-order 
correction to the triplet ground state energy is 
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FIG. 3: Occupation of the QD levels in the ground state in three 
consecutive valleys with total electron number n — 1, n, and n + 
1, respectively. A doublet- triplet-doublet spin sequence is shown. 
The variables x, (xi, X2, ■ ■ ■), and (yi, y2, ■ ■ .) denote the energy 
level spacings in the QD normalized by mean level spacing A. For 
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FIG. 4: Six distinct types of the virtual processes contribute to the 
ground state energy correction for the QD in the triplet state. Only 
tunneling processes in (or out of) the left lead are shown. In the first 
four cases spin of the dot in the virtual state S has a definite value. 
In the last two cases (e) and (f) QD spin in the virtual state has two 
allowed values: S = i and S = | with the probabilities Ws given 
by Eq. <2H . The electron structure of the virtual state corresponding 
to two allowed values of S is circled by dashed line in panel (e). 
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{similar terms for the right lead: p — » q} , (20) 



where 5 is the spin of the QD in the virtual state. One can 
easily find 5 for the first four processes, Figs.|4la)-(d), and so 
calculate the denominators for the first four terms in d20l : the 
values are |, |, |, and |, respectively. In the last two cases, 
Fig s -S e ) and (f), the QD spin in the virtual state can take two 
values, 5 = \ or |; it does so with the following probabilities 
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The corresponding contributions to the energy correction must 
be weighted accordingly. In addition, the energy difference in 
the denominators depends on 5; to account for this depen- 
dence, we introduce an additional function 



f s = 2-5(5 + 1). 
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appearing in the denominators of the fifth and sixth terms in 

Let us integrate over the continuous energy levels in the 
lead. The sum can be replaced by an integral, ^ p • • • — ► 
/ de vl{e) ■ ■ ■ , where is the density of states in the left 
lead. Taking the dot-lead contacts to be point-like, the tun- 
neling matrix elements {ifc p } depend on the momentum in 
the lead p only weakly; hence, tk p w t^L- In addition, as the 
leads are formed from 2DEG, their density of states is roughly 
independent of energy. We assume that it is constant in the en- 
ergy band of 2Ec near the Fermi surface. Then the result of 
integrating over the energy spectrum in the lead (in schematic 
form) for the first term in Eq. J20i is 
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This expression diverges, but when we calculate an observ- 
able, e.g. the shift in the position of the peak maximum 
[Eq. j!3H. we encounter the energy difference between cor- 
rections to the triplet and doublet states. The result for the 
shift is, therefore, finite: 
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In a similar fashion one can calculate the second-order cor- 
rection to the ground state energy of the doublet. The dif- 
ference of these energies at the gate voltage corresponding 
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to the peak maximum in zeroth order [Eq. dl5H . needed in 
Eq. ( 1131 . then follows. There is one resonant term, propor- 
tional to In (2Ec/T), in which the lower cutoff T appears be- 
cause of the entropic term in Eq. ill 5b . Alternatively, T would 
appear as the natural cutoff for the resonant term upon rein- 
troduction of the Fermi-Dirac distribution for the occupation 
numbers in the leads. 

For a point-like dot-lead contact, the tunneling matrix ele- 
ment is proportional to the value of the electron wave function 
in the QD at the point of contact: tk a oc ipk( r a), where a = L 
or R. Here, we neglect the fluctuations of the electron wave 
function in the large lead. Thus, the following identity is valid: 
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where the average in the denominator is taken over the sta- 
tistical ensemble. Note that by taking the ensemble average 
of both sides of d25l >. one arrives at the standard golden rule 
expression for the dimensionless conductance. 

In our calculations we take advantage of the fact that J s < 
A <C Ec and neglect terms that are of order l/£ = A/2Ec- 
Sums like 
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are split using £ fc= „ +2 • • ■ = Efc=fi+ 2 ■ • •+E fc = t +f+ 2 ' ' 
and so the last term, which is 0(l/£), is dropped. Likewise, 
expressions like 
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are of order 0(l/£), and so neglected. 

Thus, for the second-order correction to the position of the 
peak maximum, we obtain 
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where 2J S > — £| > because the total spin of the QD with n electrons is equal to 1. 

In a similar fashion one can find the shift in the position of the other peak maximum, A/n +i (l ~~ * |) • Th en ' according to 1141 . 
the difference of these two shifts yields the second-order correction to the spacing for the doublet-triplet-doublet spin sequence: 
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A potential complication is that the addition of two electrons to the quantum dot — may scramble the energy 

levels and wave functions of the QD^S^Sk Since the number of added electrons is small, we assume that the same realization 
of the Hamiltonian, Eq. (|5}, is valid in all three valleys^ 

For the second-order correction to the spacing for the doublet-singlet-doublet spin sequence we similarly obtain 
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where e^.+i — e« > 2 J s > because the total spin of the QD 
with n electrons is equal to in this case. 

Unlike the zeroth-order spacings, the second-order correc- 
tions are functions of many energy level spacings as well as 

(2) 

the wave functions at the dot-lead contact points: s = 

(2) 

U^' s (x,X,Y; {Z ka }), where x, X = (xi,x 2 , ■ ■ ■), and Y = 
(yi, ?/2, . • •) are the energy level spacings in the QD normal- 



ized by the mean level spacing A (see Fig. [5} and 

y |^fc(rq)| 2 

\\ipk(r a )\ ) 

The expressions for suggest that the main contribution 
to their fluctuation comes from the fluctuation of the energy 
level x and the wave functions {ip k (r Q )}. The other spacings, 
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X and Y, always appear within a logarithm; therefore, their Converting to dimensionless units, we find that 
contribution to the fluctuation of is small. With good 

accuracy, one can replace these levels by their mean value ,,(2) / i\ ^ ^ 1 r ^ ~n ,-,n\ 

3 F } K,s wA Z k<x}) = ^7 2^ 9 a ^ a ,s{x]{Z ka \), (33) 



Ks « !; {^» = {Z fca }) . (32) 



where 



*a,S=0 (x; {Z ka }) = (%, a + %+!,<*) 



1 3 
- ln£ + In (5 + - In a; + - ln(x - 2j) 



, " ,1 + f)- In ( 1 + ?T7 



|S( Z i-'' a + Z 9+ 1 + I '«) ln ( 1 



O (1) , (34) 



1=1 



1=1 



x + I 



®n,s=i {x; {Z ka }) = (Zj± yCC + Za+l, c 



„,,l + «)- h (l + £J 



- ln£ - - In 5 + - In 2 7 + In ( — - 1 
2 2 J \ x 



^ ( Z f-l,a + Z%+l+l,a) 
1=1 

I 



3j 



i=l 



x + z 



0(1), (35) 



where 6 — A /T. Here, the upper limit in two of the sums is 
infinity because the Fermi energy is the largest energy scale. 
In summary, the total spacing is 



where the first term is given by Eqs. J17I - J1 8i and the second 
by d33l -( l35l . The spin of the QD in the even valley, S, can 
take two values, or 1, depending on the spacing x. 



U^ a (x-{Z ka }), 



(36) 



V. ENSEMBLE-AVERAGED CORRECTION TO THE 
PEAK SPACING 

The average and rms correction to the peak spacing can now 
be found by using the known distribution of the single-particle 
quantities x and {Z ka }. In what follows, (U) denotes the 
average over the wave functions, U denotes the full average 
over both wave functions and energy levels, and P(x) is the 
distribution of the spacing x. Since (Z ka ) = 1, (& a ,s) does 
not depend on a, and the average "even" spacing is2& 



poo 

y o dxP{x)(u^ 



9L + 9R 
47T 2 £ 



dx P(x) ( $ 



2/ 



j(2) 

-a,S=0 



(37) 



(38) 



2/ 



dxP{x)(§ 



-a,S=l 



Using the asymptotic formulas 

00 

[ 



1=1 



In 1 



In 1 



x ln£, 



E 

1=1 



in 



(39) 



for £ 3> 1 in the expressions for ($ a ,s)> we find 

($«,s=o) = 2(2x-3j-l) + 2 In 5 
($ a ,s=i> = 2(-2x + 5j-l) ln£-ln^ 

valid for £, <5 ^> 1. By carrying out the integration over the 
distribution of the spacing x, the final expression is 



(40) 



C(j) = 2[-8jP (2j)+4x (2j)+5j-3], (41) 
V(j) = 3P (2j) - 1, 

where P (2j) = f™dxP(x) andx (2j) = J^dxxP(x). 
Note that Po(2j) is the probability of obtaining a singlet 
ground state while Xo(2j)/Po(2j) is the average value of x 
given that the ground state is a singlet. 

For the CI model, j = and, hence, C(0) = P(0) = 2. In this 
limit, then, the ensemble-averaged correction to the spacing is 

2E C 



(0) = 9L 



9R 



2tt 2 £ 



In- 



T 



(42) 



in agreement with previous work. 18 The magnitude here is ap- 
proximately 0.05 (gL + gji) In (2Ec/T) in units of the mean 
level spacing. 

It is convenient to relate the average change in spacing at 
non-zero J s to that at J s = 0: 



6u(j) = 



ui 2) (0) 



XC(j)+V(j) 
2(A+1) : 



(43) 
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FIG. 5: Correction to "even" peak spacing as a function of strength 
of exchange, j — J s /A, normalized by the correction at j — 0. GUE 
case with A = 1 and gL = gn- Solid: Ensemble average. Dashed: 
Ensemble average plus/minus the rms, showing the width of the dis- 
tribution. 



A 



ln£ 
In <5 



ln(2£ c /A) 
ln(A/T) 



(44) 



The dependence of 8u on j is fully determined by the param- 
eter A and the choice of the random matrix ensemble. 

Figures [5] and [6] show the results in the GUE ensemble for 
A = 1 and 3, respectively, and Fig.0shows those for the GOE 
ensemble at A — 1. In evaluating these expressions, we use the 
Wigner surmise distributions for P(x), which allow an ana- 
lytic evaluation of Po(2j) and Xo{2j). As j increases, the 
average correction to the peak spacing decreases monotoni- 
cally in all three cases. (Note, however, that our results are 
not completely trustworthy when 0.4 < j < 0.5 because in 
this regime higher spin states should be taken into account.) 
Since A depends on £ and 6 only logarithmically, the qualita- 
tive behavior of Su(j) is very robust with respect to changes 
in charging energy, mean level spacing, or temperature. 

Similarly, the dependence of 8u on A at j=jo is 



Mio, A) 



ACC?o)+X>(j ) 
2(A + 1) 



(45) 
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FIG. 6: The same quantities as in Fig.[3]plotted for A = 3. 
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FIG. 7: Ensemble-averaged correction to the "even" peak spacing 
as a function of strength of exchange, j = J s /A, normalized by the 
correction at j = for A = 1. Solid: GOE. Dotted: GUE. 



Figure [3] shows results in the GUE case for several values of 
jo- Thus, for the realistic value jo = 0.3, the CEI model gives 
an average correction to the peak spacing that is two to three 
times smaller than the CI model. 



VI. RMS OF THE CORRECTION TO PEAK SPACING 
DUE TO MESOSCOPIC FLUCTUATIONS 

Mesoscopic fluctuations of the correction to the peak spac- 
ing are characterized by the variance of . It is convenient 
to separate the average over the wave functions from that over 
the spacing x, writing 



var 



where 



dxP 



(46) 



(47) 




FIG. 8: GUE ensemble-averaged correction to the even peak spacing 
as a function of A = ln(2_Bc/A) / ln(A/T) at j = j . The curves 
correspond to jo = 0.3, 0.4, and 0.5 from the top to the bottom. 
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is the contribution due to wave function fluctuations and 

-I 2 







dxP{x)(U^f- / dxP{x)(U^ s ) 
Jo 



(48) 



is the contribution due to fluctuation of the spacing x. 

We start by considering the fluctuations of the wave func- 
tions. As the number of electrons in the dot is large, the 
distance between the left and right dot-lead contacts is large, 
r-L — vr\ 3> Xf, where Xp is the Fermi wavelength. There- 
fore, the wave functions at and yr are uncorrected, 29 



((Z kL - 1) (Z k , R - 1)} - 



(49) 



for all k and k'. The fluctuation of can then be written 
entirely in terms of the properties of a single lead: 



U 



(2) 
n,S 



_ 9l+9r/ 



(4^0' 



A 



L,S 



L,S) 



(50) 

The cross terms disappear here because, according to RMT, 
wave functions of different energy levels are uncorrelated 
even at the same point in space, 



({Z kL -l)(Z k > L -l)) = -S w 



(51) 



where /3 = 1 (j3 = 2) for the GOE (GUE) case. In fact, only the 
k = |f and k = ~ + 1 terms contribute, as one can see by using 



i 

E 

i=i 



In 2 1 - 



a; + 1 



= 0(1) 



(52) 



valid for £ 3> 1. Integrating (15 Ot over the distribution of x 
according to Eq. \A1\ (keeping in mind £, S 3> 1), we obtain 



a z(u&)) = £±& x (53) 
Z V n J /3(4tt 2 £) 2 

{41n 2 £ + [3P (2j) + 1] ln 2 5 - 4 [3P (2j) - 1] ln£ln<5}. 

In the contribution to the variance due to fluctuation of the 
level spacing x, Eq. d48i . (U-lf > s) can be taken from the previ- 
ous section. Since the average eliminates the dependence on 
the lead a, we have immediately 



(54) 



where ($l,s=o) an d (&l,s=i) for £,,S ^> 1 are given by 
Eq. (1401 . Using these expressions in Eq. d48i . we obtain 



o* = y^l* ) ( C « ln " £ + C ^ ln ' 6 + C V ^ £ In 5) . 

(55) 

Explicit expressions for the coefficients are given below once 
we reach the final result. 

The dependence on and qr of the two contributions to 
the variance is different. In particular, the contribution due to 
fluctuations of the wave functions [Eq. ( I53H is proportional to 



2 , 2 _ (_9l + gg) . (9l-9r) 2 
9l+ 9r — ^ 1 ^ 



(56) 



The first term has the same form as the contribution d55l > from 
fluctuations of x. It is convenient to write the total variance as 
a sum of symmetric and antisymmetric parts. Our final result 
for the variance is 



var 



(^>)=oJ(^ a >)+o*(Z<< 9 >) 



(57) 



where 



= (^f) 2 (%ln 2 e + <S <5<5 ln 2 5 + 5^1n£ln5) (58) 

°Z(l4 2) ) = (^f 5 ) (-Aefln 2 e + A«ln 2 <5 + ^ln£ln5) (59) 
with the coefficients {S} and {.4} given by 

<%(j) = | + 16 ( X - 1) + 64 [2jP (2j) - x (2j)] {2j [1 - P (2j)] - [1 - x (2j)] }, (60) 

SssU) = ^[3Po(2i) + l]+9P (2i)[l-Po(2j)], (61) 
S CS (j) = -1 [3Po(2j) - 1] + 24{x (2j) [1 - P (2j)] + [1 - x Q (2j)} P (2j) - 4jP (2j) [1 - P (2j)] }, (62) 

Ai(J) = I, A<5(i) = ^[3Po(2i) + l], and = [3i%(2j) - 1] . (63) 



10 



The constant \ introduced in Eq. J60> is 



VII. CONCLUSIONS 



X = I dxx P(x). 
Jo 



(64) 



For the CI model, j = and Pq(0) = xq(0) = 1; hence 

_ 4 gj + g| 



var 



j=o /3 (4tt 2 £) 2 



(ln£-ln£f 



(65) 



The first term is due to fluctuation of the wave functions at 
the dot-lead contacts, and the second term comes from the 
fluctuation of the level spacing x. The presence of the second 
term was missed in previous work [see Eq. (44b) in Ref. Hal . 
If £= 6, the first term vanishes; nonetheless, due to the second 
term the variance is always finite. 

Let us consider a realistic special case of symmetric tun- 
nel barriers, cjl = gR. 19 Then the asymmetric contribution 
vanishes, and the rms fluctuation of the correction to the 
peak spacing normalized by the average correction at j = 
[Eq. is 



zi 2) (o) 



2(A + 1) 



(66) 

Figure[9]shows this quantity plotted as a function of j for both 
GOE and GUE. Notice that (i) the rms is of the same order 
as the average, and (ii) its magnitude weakly depends on j. 
To show the magnitude of the fluctuations in the correction 
relative to its average value, we plot two additional curves in 
both Fig. |3 and Fig. [6] namely Su ± a$ u . We find that at 
the realistic value j = 0.3, the correction to the even peak 
spacing is negative for a small fraction of the quantum dots in 
the ensemble. 
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FIG. 9: The rms of the correction to the peak spacing for the sym- 
metric setup §l — gn as a function of j = J s /A normalized by the 
ensemble-averaged correction at j — 0, Eq. 1661 . The solid (dotted) 
curve corresponds to the GOE (GUE) at A = 1. The dashed curve 
corresponds to the GUE at A = 3. 



In this paper we studied corrections to the spacings between 
Coulomb blockade conductance peaks due to finite dot-lead 
tunneling couplings. We considered both GOE and GUE ran- 
dom matrix ensembles of 2D quantum dots with the electron- 
electron interactions being described by the CEI model. We 
assumed T<^A<^_ Eq- The 5 = 0, |, and 1 spin states of the 
QD were accounted for, thus limiting the applicability of our 
results to J s < 0.5A. 

The ensemble-averaged correction in even valleys is given 
in Eq. \A\l . The average correction decreases monotonically 
(always staying positive, however) as the exchange interac- 
tion constant J s increases (Figs. 15171 The behavior found 
is very robust with respect to the choice of RMT ensemble 
or change in charging energy, mean level spacing, or tem- 
perature. Our results obtained in second-order perturbation 
theory in the tunneling Hamiltonian are somewhat similar to 
the zeroth-order resultsASiii in that the exchange interaction 
reduces even-odd asymmetry of the spacings between peaks. 
While the average correction to the even spacing is positive, 
that to the odd peak spacing is negative and of equal magni- 
tude. 

The fluctuations of the correction to the spacing between 
Coulomb blockade peaks mainly come from the mesoscopic 
fluctuations of the wave functions and energy level spacing x 
in the QD. The rms fluctuation of this correction is given by 
Eqs. (I57i - d64i . It is of the same order as the average value 
of the correction (Figs. [5] and [6) and weakly depends on J s 
(Fig. [9). Therefore, for a small subset of ensemble realiza- 
tions, the correction to the peak spacing at the realistic value 
of j = 0.3 is of the opposite sign. The rms fluctuation of the 
correction for an odd valley is the same as that for an even 
one. 

We are aware of two experiments directly relevant to the 
results here. First, in the experiment by Chang and co- 
workers, 20 the corrections to the even and odd peak spacings 
due to finite dot-lead tunnel couplings were measured. It was 
found that the even (odd) peak spacing increases (decreases) 
as the tunnel couplings are increased. This is in qualitative 
agreement with the theory; see Eq. (141 i . The magnitude of 
the effect was measured at different values of the gas param- 
eter r s (and, hence J s ) as well. Unfortunately, because the 
effect is small and the experimentalists did not focus on this 
issue, one cannot from this work draw a quantitative conclu- 
sion about the behavior of the correction to the peak spacing 
as a function of J s . 

Second, in the experiment by Maurer and co-workers,— the 
fluctuations in the spacing between Coulomb blockade peaks 
were measured as a function of the dot-lead couplings with 
9l=9r- Therefore, only the symmetric part [Eq. d58l l would 
contribute to the total variance. Reference 19 reported results 
for two dots: a small one with area 0.3 /im 2 and a large one 
with area 1 /im 2 . From the area (excluding a depletion width 
of about 70 nm), we estimate that the large (small) dot con- 
tains about 500 (100) electrons. Measurements on the large 
QD found larger fluctuations upon increasing the dot-lead tun- 
nel coupling, in qualitative agreement with the theory. Though 
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the temperature was larger than the mean level spacing in the 
large QD whereas our theory is developed for T« A, the 
theory gives about the correct magnitude for the peak spacing 
fluctuations. It is inconclusive whether the data is in better 
agreement with the CI or CEI model as the fluctuations are 
roughly the same (Fig. |9j in both. In the small QD, there 
is an anomaly for the strongest coupling in the experiment - 
the fluctuations suddenly decrease. In addition, the experi- 
mental fluctuations are one order of magnitude larger than the 
theoretical estimate [Eqs. J66I and J42H . The reason for this 
discrepancy is not clear at this time. Possible contributing fac- 
tors include scrambling of the electron spectrum as the charge 
state of the dot changes, or the role of the fluctuations when 



the dot is isolated (i.e., fluctuations in U^). In order to assess 
quantitatively the role of dot-lead coupling in the Coulomb 
blockade, further experiments are needed. 
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